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Tunneling between a photo-excited p-type silicon tip and a gold surface is studied as a function 
of tip bias, tip/sample distance and light intensity. In order to extend the range of application of 
future spin injection experiments, the measurements are carried out under nitrogen gas at room 
temperature. It is found that while tunneling of valence band electrons is described by a standard 
process between the semiconductor valence band and the metal, the tunneling of photoelectrons 
obeys a Fowler-Nordheim-like process directly from the conduction band. In the latter case, the 
bias dependence of the photocurrent as a function of distance is in agreement with theoretical 
predictions which include image charge effects. Quantitative analysis of the bias dependence of the 
dark and photocurrent spectra gives reasonable values for the distance, and for the tip and metal 
work functions. For small distances image charge effects induce a vanishing of the barrier and the 
bias dependence of the photocurrent is exponential. In common with many works on field emission, 
fluctuations in the tunneling currents are observed. These are mainly attributed to changes in the 
prefactor for the tunneling photocurrent, which we suggest is caused by an electric- field-induced 
modification of the thickness of the natural oxide layer covering the tip apex. 

PACS numbers: 73.40.Gk, 72.40,+w, 79.60.-i, 72.25.Hg 



I. INTRODUCTION 

Although the photo-sensitivity of tunnel processes be- 
tween a metal and a semiconductor has been widely 
considerecUi^, the majority of these investigations con- 
cern the effect of surface photovoltage on tunneling be- 
tween a metallic tip and a planar semiconducting surface. 
On the other hand, the mechanisms governing injection of 
photoelectrons from a semiconducting tip under light ex- 
citation into a planar metallic surface are less well studied 
despite the fact that a proper understanding is important 
both from a fundamental viewpoint, and for a variety of 
potential applications. 

Photoexcited semiconducting tips can in principle be 
used as local spin injectors since circularly polarized light 
produces a spin polarized electron population via optical 
pumping 4 . The mean spin of the tunnel injected elec- 
trons can then be controlled via a change in the light he- 
licity. There exists numerous potential applications for a 
spin injector of this type, such as spin-polarized scanning 
tunneling microscopy (SPSTM)4, quantum computing 6 
and spintronics 7 -, none of which has been convincingly 
demonstrated. While it is known that a change in pump 
light hclicity modifies the tunnel junction resistance be- 
tween a ferromagnetic metallic tip and a photoexcited 
GaAs surface via the spin-dependence of the tunneling 
process^, spurious optical effects in photoexcited semi- 
conductor tips have been blamed for the apparent ob- 
servation of non-zero spin polarizations, even on non- 
magnetic surfaces 9 -! -. In order to better understand non- 
polarized photo-assisted tunneling from a semiconductor 
tip into a metallic surface a number of groups 1 ! com- 
menced by studies of the tip bias, distance and light in- 
tensity dependence of the injected current. Perhaps the 
most detailed investigations, both theoretical and exper- 
imental, were performed by the Nijmegen group 1 ^. For 



relatively small tip bias up to about 0.5 V, tunneling 
current spectra were interpreted using a model that con- 
sidered both the characteristics of the space charge layer 
formed at the tip surface, and of the tunnel barrier itself. 
While this model can account for a number of observed 
phenomena, its extension to include spin polarization ef- 
fects has met with less succes o 9 ' 13 . 

A second related area of interest is the emission 
of electrons from high aspect ratio objects such as 
semiconductin g 14 i 15 i 16 or carbon nanotube tipsii, as 
well as the influence of photo-excitation on the emitted 
current 1 ^. While the electric fields are similar to those 
used in tunneling experiments, the tip/sample (or rather 
the cathode/anode) distances are typically much larger 
(up to several /im). Under these conditions, the emit- 
ted current is usually described by a Fowler-Nordheim 
(FN) process^ 9 - from surface states to vacuum states. Un- 
less a complex sample (anode) structure is used 20 , this 
process is spin independent since the vacuum states are 
themselves unpolarized. The bias dependence of the FN 
current is of the for m 19 ' 21 

I = A(U)/( 7 U) 2 x exp(-S(V0/7V), (1) 

where -A(V) and B(V) are slowly-varying functions of 
bias. The quantity 7, larger than unity for sharp tips, 
describes the geometrical enhancement of the electric 
field 2 !. 

In the present paper we investigate the mechanisms 
which govern photo-electron injection from a semicon- 
ducting tip. The tip bias value, as large as -3 V, is suf- 
ficient to induce an increase in the injected current by 
FN-tunneling, but not large enough to enable observa- 
tion of current oscillations related to quantized states in 
the tunneling gap2£. Here, we are not interested in spin- 
polarized injection so silicon tips and gold surfaces are 
used. With respect to a number of previous works on 
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FIG. 1: Injection geometry: a semiconducting tip, placed at 
the end of an AFM-like cantilever, is excited by light imping- 
ing on its rear planar face. 



tunnel injection of photocarrier o 9 ' 10 i 12 ' 13 i 23 , the present 
situation is simpler since: i) the light excitation is inci- 
dent on the rear planar surface of the tip, which itself is 
situated at the end of a very stiff cantilever (see Fig. []}. 
Consequently, the photoelectrons diffuse from the rear 
of the tip (where they are created) to the tip apex, and 
the number of photoelectrons reaching the tip apex is in- 
dependent of tip bias. This is at variance with studies 
where the light excitation is incident on the front face 
of the tip ape x 12 ' 13 , and the injected photocurrent de- 
pends on tip bias via the bias dependent depletion layer 
width 2 !. The experiment reported here thus resembles 
an atomic force microscope (AFM), with the exception 
that tip/sample distance is controlled using the tunnel 
current 25 , and ii) the use of p + doped tips guaranties 
that the conductivity of the tip is high and that most 
of the applied bias is dropped across the tunnel barrier 
itself. The experimental results can then be understood 
using a relatively simple theory, iii) the tip surface bar- 
rier (Vf,) can be neglected for V > V&, and iv) the ex- 
periments are performed in an inert gas atmosphere. Al- 
though this has intrinsic drawbacks compared with ultra- 
high vacuum conditions, in particular related to possible 
tip pollution and current instabilities caused by changes 
in the native oxide layer thickness, the choice is deliber- 
ately made in order to simplify and extend the scope of 
application of future SPSTM studies^. 

The dependence of the tunneling dark current and pho- 
tocurrent as a function of bias, tip/sample distance and 
excitation light intensity is investigated. It is found that 
i) the dark current obeys a standard tunneling process 
between the valence band of the semiconductor and the 
metallic local density-of-states, and ii) the tunneling pho- 
tocurrent occurs directly from the conduction band of 
the semiconductor to the metallic states through a FN- 
like process and not, as found by the Nijmegen group 
for GaAs tips^, via a standard tunneling process from 
midgap states at the semiconductor surface. The impli- 
cations of these findings for spin-polarized tunneling are 
discussed. 
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FIG. 2: (a) Schema of the band structure of the semiconduc- 
tor (left) and of the metal (right) for a tip bias V. Arrow 
1 represents the dark current arising from a tunneling pro- 
cess between the semiconductor valence band and the metal. 
Arrow 3 represents the current arising from a tunneling pro- 
cess between localized mid-gap semiconductor states and the 
metal, as proposed in the Nijmegen modeU^. Arrow 2 shows 
the observed FN-like tunneling process of photoelectrons from 
the semiconductor conduction band, (b) Tunnel barrier shape 
with (solid line) and without (dotted line) image charge ef- 
fects. 



II. THEORETICAL CONSIDERATIONS 

The following section describes tunnel injection be- 
tween a negatively biased, p + doped photoexcited semi- 
conducting tip and a metallic surface. Shown in Fig. 
[5^ is a schema of the band structure of the tip and 
metallic surface at applied tip bias, V. Light excita- 
tion from the rear of the tip results in the formation of 
a steady state population of photoelectrons in the con- 
duction band. Since the tip is highly conductive, the 
semiconductor band structure is independent of bias and 
is simply shifted with respect to the metallic one by qV, 
where q < is the electronic charge. 

Shown in Fig. is a schematic of the shape of the 
tunneling barrier. Because of image charge effects^ 7 -, the 
tunneling barrier at the energy of the tunneling electrons 
extends between d\ and d-i-, and its width Ad — (I2 — d\ 
is smaller than the physical tip/sample distance, d. 

Because the applied tip bias is negative, the dark cur- 
rent and the current of minority photocarriers are de- 
scribed by a tunneling of electrons from the semiconduc- 
tor to the metal. The incremental tunnel current from 
an occupied energy level E in the tip to an empty level 
in the metal at the same energy is given b y 27 ' 28 



Sit = K S p s p rn exp(— nAd) , 



(2) 



where p s and p m are the density of states at E of the 
semiconductor and of the metal respectively, and S is 
the tip surface area. The quantity K is proportional to 
the tunneling matrix element, and depends on the elec- 
tronic orbitals of tip and surface via which tunneling oc- 
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curs. The inverse distance for the exponential decay of 
the tunneling current is given by2£ 



(3) 



where m is the mass of the electron in vacuum. Here, the 
zero of energy is taken at the Fermi level of the semicon- 
ductor. In the work of Simmons 27 , E is the sum of the 
potential energy of the tunneling electrons and of their 
kinetic energy oriented perpendicular to the surface E z . 
In the present case, since tunneling from a tip is essen- 
tially a unidirectional process 2 ^, the tunneling electrons 
have a velocity along the tip axis and E z — E. Eq. [^ex- 
presses the highly simplifying result that tunneling across 
a barrier of complex spatial shape mostly depends on the 
spatial average of the barrier, defined as 



1 

Ad.,,, 



(z)dz. 



(4) 



A. Tunneling photocurrent 

In general the tunneling current obtained from Eq. [2l 
Eq. [3] and Eq. [4] is given by the expression 



Iph = A P h exp(-K ph d) 



(5) 



The prefactor A p h is found by replacing p s by the con- 
centration N e of photoelectrons (a quantity proportional 
to the excitation light power): 



Aph = KphSNePn 



(6) 



For a standard tunneling current in the absence of image 
charge effects one has 4> = (f) a —qV*/2. Since the photoex- 
citation is only slightly above gap E is the semiconductor 
bandgap energy E g and K p h is given by 



Kph 



(7) 



where \o = </"o — E g is the affinity of the tip surface. The 
effective tip bias V* is given by 

V* = V + V = V + -(4> -<t> m ) (8) 
q 



between tip and 
6 m ~ —0.6 eV for 



since there is an electric field at V 
surface if <p m is different from (f>o (4>_ 
clean surfaces of silicon and gold) . 

At larger tip bias, tunneling from the semiconductor to 
the metal can occur via vacuum states (see Fig. [2^,, arrow 
2). This process resembles FN emission with three differ- 
ences with respect to Eq. [T]and will therefore be denom- 
inated FN-like. Firstly, in cases where the tip/sample 
distance is small compared with the radius of curvature 
of the tip, no field enhancement occurs at the semicon- 
ducting tip (i.e. 7=1 in Eq. [T]as explained in Appendix 



1). A second difference with the field emission process 
is the replacement of V by V* given by Eq. [5] (in field 
emission studies 0o — 4>m is negligible compared to the 
applied tip (cathode) bias and can be ignored). Finally, 
since the tunneling photoelectrons are distributed over 
a narrow energy range at the bottom of the conduction 
band, it is not necessary to perform an integration over 
energy. Unlike Eq. [TJ the prefactor is independent of tip 
bias so that Eq. [5] is still valid, albeit with a modified 
value of Kph ■ 



FN-like tunneling without image charge effects 



In order to observe FN-like tunneling, qV* must be 



larger than 



Eg, so that the threshold bias Vth is 



qV th = 4> m - E g 



(9) 



The spatial average 4> of the barrier potential 3 -^ is then 
equal to (4>o + Eg) /2. The tunneling photocurrent is ob- 
tained from a calculation of the quantum tunneling prob- 
ability at E gi and is given by Eq. [5] with 



Kph 



3/2 

m Xq 
h qV* ' 



(10) 



If Vq is negligible with respect to V , Eq. [10]indicates that 
a plot of Iph in logarithmic units as a function of l/V 
should be linear. The slope of this line should decrease 
when the d is decreased. For a bias lower than Vth, the 
photocurrent has a smaller value obtained from Eq. [5] 
and Eq. [7] 



2. FN-like tunneling with image charge effects 

Inclusion of image charge effects requires a modifica- 
tion of Eq. [TO] It is not a bad approximation to ap- 
ply the treatment of Simmons for tunneling between two 
metals 27 , since with respect to the static dielectric con- 
stant of vacuum (e = 1), that of silicon (e = 13) can be 
approximated by that of a metal (e — > oo). The image 
charge potential is given by 2 ^ 



(z) = 4> - 



qV*z 



1.15 



\d 2 



z(d — z) ' 



(11) 



where A describes the magnitude of the image charge 
effects: 



A = 



q 2 ln(2) 
Sireead 



(12) 



Here eo is the permittivity of free space. For qV* suf- 
ficiently large 31 and d\ <§; d, one finds that the inverse 
distance for tunneling is now given by 

3/2 



m X (1 ~ 5.6A/xo) 
qV* 



v/l + 5.6A/ X o 



aqV*/xo, 
(13) 
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where 



Xo 



1.2 



2.3 



■In 



l-5.6A/xo \d 1 (d-d 2 ) 



did 



(14) 



is weakly dependent on V* via the logarithmic term. An 
approximate^ but physically meaningful expression for 
the threshold bias Vth for FN-like tunneling is found by 
setting di = d which gives 



qVth = 4>n 



E a - 5.6A, 



(15) 



which decreases with decreasing d unlike the value given 
by Eq. [9] As expected, if A <C Xo, Eq. [TBI reduces to Eq. 
[TU]and Eq. [15] reduces to Eq. 

The inclusion of image charge effects modifies the pho- 
tocurrent in two ways: i) at small bias: because of the 
dependence of n p h on d in Eq. [131 the slope of log(/ p / l ) 
as a function of 1/V is no longer proportional to d un- 
less this distance is large. The effective tip affinity Xo 
is now given by Xo 3 = (Xo + 5.6A) (xo — 5.6A) 2 , and ii) 
if aqV* « xo + 5.6A the lowering of the barrier due to 
image charges causes an excess photocurrent. 
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FIG. 3: Schema of the atomic force microscope-like experi- 
mental setup. The sample is fixed on a piezoelectric tube, 
enabling scans in the plane of the sample. The cantilever 
holding the tip is fixed on a second piezoelectric tube enabling 
motions in the vertical direction. 



III. EXPERIMENTAL DETAILS 



B. Tunneling dark current 

As represented in Fig. [2] the dark current is pre- 
sumed to occur from tunneling processes between oc- 
cupied states in the semiconductor valence band and 
empty states in the metal at the same energy. At va- 
lence band energies, image charge effects^ and FN-like 
tunneling can be neglected. An energy-independent den- 
sity of states for the metal is assumed^, and the usual 
form of p s for a bulk semiconductor is used, p s {E) = 
l r2m^i3/2 v /__g w here h is Planck's constant, m* is 
the effective electronic mass, and E < is the energy of 
the tunneling electrons in the valence band. The value 
of the dark current is obtained by integrating Eq. [5] be- 
tween the Fermi level of the semiconductor and that of 
the metal, 

hark = A dark / \ r ^E exp(-K dark d)dE (16) 

J-qV 

where Kdark and Adark are given respectively by 



2m 

K dark =2\j—{4> -E-qV*/2), (17) 



and 



Adark — E-dark 



Sp m ,2m* 



27T2 L h 2 



1 3/2 



(18) 



The dark current is the sum of several exponential con- 
tributions and generally does not have a simple exponen- 
tial dependence as a function of distance. 



Experiments were conducted in a tailor-made system 
enabling STM as well as AFM investigations at room 
temperature, and consistent with the geometry of pho- 
toelectron injection described in Fig. Q] This sys- 
tem, shown in Fig. [3] operates in an electromagneti- 
cally shielded inert gas environment. Worth mentioning 
are the following particularities: i) the excitation laser, 
(wavelength 780 nm, power 10 mW) is focused on the 
back planar face of the cantilever, yielding a spot of di- 
ameter 50 /im opposite the tip position^ 4 , ii) a beam- 
splitter situated between the laser and the cantilever 
sends the light reflected from the cantilever surface onto a 
quadrant photodiode, thus enabling simultaneous atomic 
force, dark current and photocurrent measurements, iii) 
two piezoelectric tubes are used. Scans in the plane of 
the sample for imaging purposes are performed using the 
four-electrode PZT tube on which the sample is mounted, 
while the tip/sample distance is changed by a second two- 
electrode PZT tube on which the block holding the can- 
tilever and tip is mounted. This ensures that the light 
spot is stationary with respect to the location of the tip, 
iv) for current measurements the sample is grounded and 
the bias is applied to the tip. The tip/sample current is 
monitored using a high gain/low noise amplification cir- 
cuit placed as close as possible to the tip. After installing 
the sample and cantilever in an air ambient, a metallic 
grounded hood is placed over the experiment and dry ni- 
trogen is blown for 20 minutes prior to, and during, the 
measurements . 

Commercially available tapping-mode silicon 
cantilevers— were used in all experiments. The 
use of relatively stiff cantilevers (of nominal force con- 
stant 16 N/m) minimizes the effect of atomic forces on 
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FIG. 4: Measurement sequence (see text). 
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tip/sample distance, and thereby improves the control of 
this distance using the tunnel current in a feedback loop. 
Silicon tips were situated at the free end of the cantilever 
and ranged in height between 5-7 fj,m, and were p-type 
doped in the range 10 20 cm" 3 . The samples consist of 
polycrystalline gold films, of thickness approximately 
60 monolayers, deposited on (111) silicon substrates by 
electrochemical methods^. 

Care was taken in order to minimize fluctuations in 
the tunnel current due to (photo)electro-chemical reac- 
tions which can induce changes in the tip characteristics. 
Anodic tip potentials (forward bias for p-type tips) were 
avoided since they are known to rapidly generate a thick 
layer of oxide which prevents the observation of tunnel 
currents. Several tip passivation treatments were tried in 
order to render the tip surface inert, including removal 
of the native oxide with hydrofluoric acid. However, the 
native oxide layer yielded the most inert surface, and the 
results reported here were all obtained with naturally ox- 
idized tips. As will be seen, the presence of residual water 
can still induce electrochemical changes in the tip surface 
which produce instabilities in the tunneling current. 

The measurement sequence is shown schematically in 
Fig. During T\ w 0.5 s (i.e. much longer than the time 
constant of the feedback loop) the pump laser is switched 
off and no data acquisition takes place. The tip/sample 
distance is stabilized at V set = — 2 V by imposing a dark 
tunneling current Idark — Iset using a standard feedback 
loop. Spectroscopic measurements of the tunneling cur- 
rent as a function of tip bias are then performed in a 
rapid sequence of two acquisitions over periods T 2 and 
T3, each of which lasts 25 ms. Since these scans induce a 
change in current, the feedback loop is opened for T 2 + 
T3 = 50 ms, a period significantly shorter than time con- 
stants over which drift in d occurs. During T 2 the laser is 
switched on, while during T 3 it is switched off. The tun- 
neling dark current spectrum is that obtained during T3, 
while the tunneling photocurrent spectrum is that ob- 
tained by taking the difference of the spectra measured 
during T 2 and T 3 . Most of the results presented in this 
paper were obtained after only one such acquisition (i.e. 
without averaging over several spectra). 

Using this procedure, the dark current, the photocur- 



FIG. 5: Dependence of the atomic force between tip and sam- 
ple as a function of the dark current value imposed during T\ . 
For each value of Idark ten successive measurements were ob- 
tained, each of which is shown in the figure. The arrows, a-e 
correspond to the curves shown in Fig. [8] The abrupt change 
near Idark = 2 nA is due to mechanical tip/sample contact. 

rent, and the atomic force between tip and surface (mea- 
sured using the quadrant photodiode), were simultane- 
ously monitored for increasing tunneling dark currents 
during Tj., corresponding to decreasing tip/sample dis- 
tances. In order to check reproducibility and possible 
current instabilities, each scan was performed ten times 
in identical conditions. 



IV. RESULTS 

Fig. [5] shows the variation of the atomic force between 
the tip (shown in the inset of Fig. [5]) and the sample, 
as a function of the tunneling dark current set during 
time T\. Shown by arrows in the figure are the selected 
dark current values for which tunneling spectra will be 
presented below. The atomic force stays approximately 
equal to zero up to a dark current of the order of 2 nA, 
above which it abruptly becomes repulsive because of 
mechanical contact with the surface. This abruptness 
as well as the fact that a tunneling current is observed 
before the onset of a non-zero atomic force are due to 
the relatively large stiffness of the cantilever. Points a, 
b, and c correspond to the out of contact regime, while 
point d illustrates the near-contact regime. At point e, 
there is no doubt that gold indentation occurs since this 
has been observed 37 for forces as small as several nN. 



A. Instabilities 

Shown in Fig. [5^ are ten dark current spectra at point 
c in Fig. [5] The corresponding tunneling photocurrent 
spectra are shown in Fig. [5Jd. Both currents strongly in- 
crease with bias and no surface photovoltage is observed 
(no photocurrent is detectable below 1.15 V). 
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FIG. 6: Ten successive measurements of (a) dark current and 
(b) photocurrent as a function of tip bias for fixed tip/sample 
distance. In (a), three curves (designated by arrows) strongly 
differ from the other curves at V se * and were not considered in 
the analysis, (c) The linear dependence of the photocurrent 
of the remaining seven curves on logarithmic plot of the pho- 
tocurrent versus 1/V reveals a FN-like tunnel process. Shown 
(inset) is a scanning electron microscope image of the sharp 
tip. 



The differences in the curves reveal instabilities in the 
tunneling process. These instabilities are larger for the 
tunnel photocurrent, which is consistent with the fact 
that immediately before the measurement, the feedback 
imposes a constant dark current Idark at V = V se t ■ They 
manifest themselves as abrupt temporal changes between 
a relatively small number of well-defined values which 
can differ by up to a factor of three. Such behaviour has 
been reported in works on field emission 1 ^, on transport 
in metal-oxide-silicon structure o 38 i 39 , and in tunneling 40 . 
Since the delay between the instabilities is of the order 
of a fraction of a second (i.e. much larger than T2 + T3), 
no apparent jump is seen within a given curve. 

Of the ten dark current spectra in Fig. [S^, seven curves 
approximately coincide. The other three curves are dis- 
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FIG. 7: Blunt tip: Ten successive measurements of (a) dark 
current and (b) photocurrent as a function of tip bias for 
fixed, large tip/sample distance. Shown (inset) is a scanning 
electron microscope image of the tip. 



carded since the current value at V se t— -2 V differs from 
the setpoint value. As seen from Eqs. [5]fTUl Eq. [THl and 
Eq. 1181 instabilities can originate from fluctuations in d, 
(c/>o) or ((f> m ), or in the prefactors A p h and Adark- Insta- 
bilities in the prefactors can be caused by fluctuations of 
the tip surface S or of the tunneling probabilities K p h or 
Kdark (it is assumed that N e and p m are stable). 

Instabilities in d can be ruled out based on results ob- 
tained with a blunt tip made by mechanically removing 
the tip apex. The tip is almost ideally flat, with an end 
diameter of f=s 2 /im (see inset of Fig. [7]). For a blunt tip, 
tunnel processes are averaged out over a large area and 
should therefore be relatively insensitive to local changes 
in geometry and tip/sample surface composition. The 
spectra shown in Fig. [7] were obtained by fixing the tun- 
neling dark current during T\ at the same value used for 
the sharp tip. The dark current spectra are essentially 
indistinguishable, while differences between photocurrent 
spectra are comparable to the noise. The same qualita- 
tive results are also obtained for all tip/sample distances. 
This proves that any tip/sample distance noise or drift is 
effectively compensated by the feedback loop, and there- 
fore that fluctuations in d are not the root cause of cur- 
rent instabilities measured with the sharp tip. 

Fig. [5b shows I p h scans plotted in logarithmic units 
as a function of 1/V for the sharp tip. While I p h varies 
by up to a factor of 3, the slope (which is related to 0o 
according to Eq. fTTI)) of the curves is the same within 
experimental error. Variations in c/>q can therefore be 
ruled out. Fluctuations of Adark can also be excluded 
using the following reasoning: a change in Adark will in- 
duce a change in d via the feedback loop, and thus a 
change in the shape of I dar k(V)^. This change is not 
clearly observed (see Fig. |6^i). The only remaining pos- 
sible sources of current fluctuations are changes in A p h, 
and to some extent <p m . Since Adark (and thus S) are 
constant, changes in A pn must originate from changes of 
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the prefactor K p h of Eq. [TUJ 

The same analysis, performed for all distances, shows 
that: i) at large distances (point a of Fig. [S]) the instabil- 
ities are relatively small. This may be because the local 
electric field is smaller and thus (photo)electro-chemical 
modification of the surface oxide layer is less important, 
ii) for reduced distances before mechanical contact, the 
instabilities mostly concern K p h, iii) for measurements 
performed under mechanical contact (point e of Fig. [5]) 
the instabilities are very large. This is again consistent 
with the notion that greater electric fields result in an 
increase in oxide modifying (photo)electro-chemical pro- 
cesses. These findings are confirmed by the quantitative 
analysis presented in section V.B, where a possible mi- 
croscopic interpretation will be proposed. 



B. Tunneling processes 

The bias dependence of the dark current for different 
distances is shown in Fig. [St, while that of the tun- 
neling photocurrent is shown in Figs. [5t> and [St. For 
clarity, only curves for which the photocurrent is largest 
arc presented, although the results are typical of a much 
larger number of experiments. It is apparent that the 
tunneling dark current and the tunneling photocurrent 
originate from distinct mechanisms. The tunneling dark 
current tends to present a flatter spectral response than 
the photocurrent when plotted against tip bias on a log- 
linear scale. Furthermore, the slope of the tunneling dark 
current is relatively insensitive to a change in distance 
whereas the slope of the tunneling photocurrent spectra 
increases greatly with an increase in distance. 

For curves a, b, and c in Fig. [St, the dependence of 
ln(Iph) as a function of 1/V is quasi-linear over two or- 
ders of magnitude, consistent with Eq. [10] and/or Eq. 
1131 in the small bias regime. The tunneling photocur- 
rent behaviour is thus consistent with the FN-like process 
described in section II whereas it will be seen that the 
tunneling dark current is better described by a standard 
tunneling process from the semiconductor valence band. 

Qualitatively, the behavior of the threshold bias for 
the onset of the FN-like photocurrent, Vth, as a function 
of distance can only be understood by including image 
charge effects. Experimentally, V t h is estimated at -2.6 
V, -1.6 V, and -1.15 V for curves a, b, and c in Fig. [5b 
respectively. Eq. [9] which neglects image charge effects, 
predicts Vth ~ —4.3V (using <p m — 5.4 eV for a clean 
gold surface) which is too large and independent of bias. 
On the other hand, the experimentally observed thresh- 
olds can be understood by using Eq. [15] with tip/sample 
distances of 1.56 nm, 1.0 nm and 0.9 nm for curves a, b 
and c respectively. 

For curve d of Fig. [St, \n{I p h) above 1.5 V is larger 
than the linear extrapolation from lower bias (dotted 
line). Despite being in the near contact regime (see point 
d of Fig. O , the excess current is not caused by electro- 
static forces acting on the cantilever, as seen from the 
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1 /tip bias, |V|" 1 (V 1 ) 



FIG. 8: Plot of the dark current (a) and of the photocurrent 
(b) as a function of V for the sharp tip shown in Fig. [S] (c) 
Photocurrent spectra versus 1/V show exponential behaviour 
over almost two orders or magnitude. This, and the increase 
in slope with increasing tip/sample distance (from curve d to 
curve a), indicate a FN-like tunneling process for photoelec- 
trons. The excess photocurrent observed at high tip bias on 
curve d is also consistent with this interpretation. Solid black 
curves are fits to the data using Eq. [jljfor the dark current 
and Eq. [13] for the photocurrent (see section V). 



absence of a bias-dependent atomic force. Rather, it is 
the result of a lowering of the tunnel barrier due to image 
charge effects according to Eq. [T3] 

In the case of curve e of Fig. [Sb and Fig. [St, the pho- 
tocurrent is large. It is estimated that several 10 -3 of the 
total number of photoelectrons created in the tip are in- 
jected into the metal. No linear dependence of ln(/ p ^) as 
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a function of 1/V is observed over the whole bias range. 
In Fig. [5b, curve e above a tip bias of 0.5 V shows a sim- 
ple exponential dependence on tip bias, consistent with 
the vanishing of the tunnel barrier—, and will be dis- 
cussed further below. Since the tip is now in mechanical 
contact with the surface (see Fig. [5]), the width of the 
tunneling gap is constant and equal to the oxide thickness 
and only the area of the surface indentation, S, changes 
with set point. When in mechanical contact, it is in- 
deed observed experimentally that I p h{V) and Idark{V) 
change only by a multiplicative factor with a change in 
set point^ 2 -. 

The dark current bias dependence is found to be ap- 
proximately exponential at large tip/sample distances 
(see curves a, b and c of Fig. [Sh.). This dependence, 
which occurs over a small bias range, can be understood 
by expanding Eq. [17] to first order in V . 



TABLE I: Fitting parameters for the curve fits shown in Fig. 
[51 and for six spectra measured at point c in Fig. [S] A p h 
influences only the photocurrent, while d, <f>o, and <f> m are 
shared between the two according to Eq. Eq. [T3]and Eq. 
1161 Adark is an indirectly shared parameter via the feedback 
loop. 



spectrum A dark (A) A ph (A) d (nm) <p 4>, 



a 


656 


1.9 x KT 6 


1.52 


4.4 


4.5 


b 


171 


4.1 x KT 6 


1.29 


5.1 


6.7 


c(i) 


4.4 


6.5 x 10" 9 


1.27 


4.2 


5.3 


c(ii) 


5.4 


7.8 x 1(T 10 


1.18 


3.8 


3.3 


c(iii) 


5.4 


5.6 x 10" 9 


1.23 


4.1 


4.7 


c(iv) 


5.1 


3.6 x 10" 9 


1.25 


4.1 


4.9 


c(v) 


3.9 


2 x KT 10 


1.17 


3.7 


3.2 


c(vi) 


4.5 


4.5 x 10" 9 


1.21 


4.2 


4.9 


d 


1 


1.85 x KT 8 


1.05 


4.6 


5.7 



INTERPRETATION 



Quantitative analysis of the I p h (V) and /dark (V) 
bias dependences 



A simultaneous fitting procedure, with equal weight, 
of I P h{V) and of Idark{V) was used that consisted 



of selecting initial values 
A,, 



of 



^0, <Pm, Adark 

iph, and in finding the set of values for these quanti- 

r i 2 

ties which minimizes R = Yl Iph (Vi) — lffl c ' (Vi) + 
J2 [hark (Vi) - I^ l r c k {Vi)] where the index i labels the 

i 

individual data points. The calculated values (Vi) 

and I^rk (^») were obtained using Eq. [5] and Eq. [TBI for 
the photocurrent, and Eq. [TC] for the dark current (for 
which the integral was evaluated numerically) . Although 
there are five fitting parameters, the fit is severely con- 
strained by the fact that both currents depend on d, 4>o, 
and 4> m explicitely, as well as indirectly on Adark via the 
feedback process. This co-dependence allows for an inde- 
pendent determination of d and 4>q which is not possible 
in many works on field emission 1 -. The fits, shown as 
solid lines in Fig. [5J are in satisfactory agreement with 
both the experimental results and the qualitative analy- 
sis of the preceding section. The excess photocurrent in 
curve d at large bias is also accounted for. For the dark 
current, the small but systematic discrepancy near the 
threshold may be attributed to the surface barrier V& at 
the tip apex. If \ V\ < V , the tunneling dark current is re- 
duced because majority carriers must also tunnel across 
the space charge layer, which is estimated^ to have a 
width of order 20 nm. 

The values of the parameters obtained using the above 
procedure, and corresponding to the various curves of 
Fig. [SJ: and Fig. [SJare summarized in Table 1. 



and 



B. Origin of current instabilities 

Examination of the results of Table 1 concerning curve 
c of Fig. [6b confirms the qualitative analysis of section 
IV. A and suggests the following comments on the nature 
of the current instabilities: 

i) d takes relatively stable values around 1.2 nm. The 
magnitude of A is of the order of 0.4 eV. According to 
Eq. [T5] the lowering of the barrier due to image charge 
effects cannot be neglected. 

ii) The values obtained for (f>o and <j) m are reasonable: 
4>q lies between 3.7 eV and 4.2 eV and <p m lies between 
3.3 eV and 5.2 eV. 

iii) The largest instability concerns the quantity A p h, 
which fluctuates by a factor of 30 between the various 
curves. m also fluctuates by about 1.5 eV. In con- 
trast, the fluctuations of Adark, 4>o an d d are significantly 
smaller. Since Adark (and thus S) does not fluctuate, 
variations in A p h are related to changes in K p h- 

The instabilities are presumably related to the fact 
that the experiments take place in a slightly humid, 
gaseous environment. With a clean tip and sample sur- 
face in vacuum the only possible mechanism for instabil- 
ities is the desorption of gold or silicon atoms under the 
effect of the electric current. This process may explain 
the observed changes in <j) m if desorbed silicon atoms are 
adsorbed onto the gold surface, but cannot explain the in- 
stabilities of K p h which themselves are most reasonably 
associated with the presence of an oxide layer covering 
the tip apex. The thickness (d ox ) of the oxide is not 
known, but it is assumed to be smaller than the smallest 
value found for d in Table I « 1 nm (the thickness^ of 
the planar natural oxide ranges between 0.5 and 1.5 nm). 

Since esi0 2 = 4 the electric field is mainly dropped 
in the tunneling gap between the oxide surface and the 
metal. The relevant distance for dark current tunneling 
is still the total distance, d, between the silicon surface 
of the tip and the metal, whereas for FN-like tunneling 



9 



which is an electric field driven effect, it is d—d ox . In the 
simplest case where image charge effects are neglected, 
this results only in a modification of A p h'. 



Aph ~ Aph SXp^Kphdo 



(19) 



where n p h is defined in Eq. [101 Eq. QjUprovides a possible 
explanation for the observed variations in Aph as any 
small variation in d ox translates into large fluctuations 
in the photocurrent prefactor, A* h . In the case of the 
blunt tip, variations in A p h are averaged out over a large 
surface area: when I p h is plotted against 1/V for different 
tip/sample distances as shown in Fig. the extrapolated 
curves pass through the same point at 1/V = indicating 
that Aph is nearly identical for all distances. It is partly in 
consequence of this that fluctuations in the photocurrent 
(see Fig. [Tfc) are minimal. 

The presence of an oxide layer and the inclusion of 
image charge effects may alter the average 4> over the 
oxide and vacuum barriers and therefore n p h, but a de- 
tailed analysis is unreasonable given the additional ad- 
justable parameters (e.g. oxide thickness, dielectric con- 
stant and effective mass) which are not known if the oxide 
is ultra-thin. The impact of an oxide layer on the results 
was evaluated using a perturbation-like approach which 
consisted of introducing an extremely thin oxide (0.1 to 
0.2 nm) into the equations of gradually increasing thick- 
ness, and repeating the numerical fit described above. 
While Adark was increased significantly, the most phys- 
ically meaningful parameters (d, </>o) were unchanged so 
the main conclusions of the present work are unaffected 
if, as assumed above, d ox < 0.5 nm. 

Investigations of instabilities in field emission^ sug- 
gest a possible microscopic mechanism for fluctuations 
in d ox : the tip is thought to be covered by a layer of 
adsorbed molecules or ions (water, oxygen, etc..) that 
diffuse into (and out of) the oxide layer under the effect 
of an electric field, thereby changing its effective thick- 
ness. In view of the fact that tunneling occurs via only a 
handful of atoms/molecules at the tip apex, the observed 
sharp changes in tunneling current between certain well- 
defined values can be attributed to the adsorption or des- 
orption of individual molecules. Other explanations in- 
volving charge build-up at local defects in the oxid o 38 i 40 
can be discounted since this would imply a change in ef- 
fective barrier height which, in the majority of cases, is 
not observed. 



C. Effect of tip/sample distance 

From Table U d decreases as expected from curve a to 
curve d. The change in slope of the curves is greater 
than an exp(— d/V) dependence and arises from the de- 
pendence of A and thus K p h on distance. 

Adark decreases steadily by almost 3 orders of magni- 
tude when d is decreased. This might be attributed to a 
progressive reduction in S with decreasing distance, but 
a similar trend in A p h which also depends linearly on S, 
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FIG. 9: Plot of the photocurrent as a function of 1/V for the 
blunt tip shown in Fig. [7] The shared origin indicates that 
Aph is nearly identical for all four curves taken at different 
tip/sample distances. 



is not clearly observed perhaps due to the fluctuations in 
Aph- It may be that Adark is under-estimated at small d 
due to the presence of the oxide layer which is neglected 
in the model (see the preceding sub-section). At small 
d the relative effect of the oxide is larger and will result 
in an artificial reduction of Adark- Further studies are 
necessary to clarify whether this is the case, or whether 
the trend in Adark is due to a change in Kdark (i-e. the 
orbital overlaps) and/or S. 

In the case of curve e of Fig. [SJa, taken under mechan- 
ical tip/sample contact, the exponential bias dependence 
of the tunneling photocurrent can be explained assuming 
a complete vanishing of the tunnel barrier. Assuming 
that d is smaller than half the de Broglie wavelength of 
the electron, and taking a square barrier, it has been rec- 
ognized (see Fig. 2.6 of Ref. |H) that the probability 
of ballistic transmission is an exponential function of the 
barrier height. Taking the barrier height to be the spa- 
tial average calculated using Eq. [4] the probability for 
ballistic transmission is of the form 



T = exp {—(</>— E g )/(f)*} , 



(20) 



where <jf is some characteristic energy. Again using Sim- 
mons' approximation for the image charge potential, Eq. 
ITTT the ballistic photocurrent is of the form 

I ph = A ph exp{-( Xo /2 + 2.8\-aqV*/2) /</>*}. (21) 

The barrier given by Eq. [IT] only vanishes if A is larger 
than about 0.5 eV (this assumes for simplicity that </>o = 
<f> m = 4 eV). Taking an oxide thickness of 0.5 nm as 
suggested by section V.B, and esio 2 = 4, Eq. [121 yields 
A = 0.25 eV which is too small. However the dielectric 
constant of the ultra-thin oxide layer may be smaller than 
the bulk value, so A » 0.5 eV is reasonable. Using this 
value and taking the oxide effective mass as given in Ref. 
[451 a de Broglie wavelength for V = -2 V of the order 
of 3 nm is obtained, which is indeed larger than twice 
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FIG. 10: Dependence of the photocurrent on excitation light 
power. For the three values of the tip bias, this dependence 
is linear. 



d ox . Similarly, using Fig. 2.6 of Ref. |2g, (j>* is calculated 
to be 0.025 eV. In view of uncertainties in the value of 
d ox and in the shape of the barrier, the agreement with 
the experimental value of 0.1 eV can be considered as 
satisfactory. 



D. Comparison with other works 



In the Nijmegen model, all tunneling occurs via midgap 
surface states uniformly distributed in energy and is for- 
mulated for low tip biasi 2 -. This mechanism, illustrated 
in Fig. [2] (arrow 3), can be excluded here for two rea- 
sons: i) the most significant tip bias dependence of the 
tunneling current in the Nijmegen model arises via the 
bias dependence of the surface depletion layer width 24 , 
and is too weak to explain the strong variation observed 
in this work. Extension of the Nijmegen model to include 
the bias dependence of n p h , given by Eq. is also unable 
to explain the results, and ii) the tunneling photocurrent 
in the Nijmegen model depends logarithmically on the in- 
cident light power in contrast with the linear dependence 
observed here (see Fig. [TU]) . 

Besides the tip bias values used in Nijmegen work, 
which are too small to observe FN-like tunneling, their 
use of GaAs tips is also significant. In defect free GaAs 
the surface recombination velocity is orders of magnitude 
larger than that for defect free silicon 4 ^. Consequently, 
the characteristic time for trapping at midgap states may 
be faster than the characteristic tunneling time from the 
conduction band so tunneling will proceed via midgap 
states. The opposite may be true in silicon so that FN- 
like tunneling is favored. 



VI. CONCLUSION 

Photo-assisted tunneling between a p-type silicon tip 
illuminated from the rear and a metallic gold surface has 
been studied for absolute tip bias up to 3 V. The addi- 
tional tunneling current induced by light excitation, due 
to injection of photocarriers from the tip into the sur- 
face, is large and comparable with the dark current. Its 
magnitude depends linearly on the light excitation power 
and no surface photovoltage is observed. 

The tunneling dark and photocurrent spectra are dis- 
tinct and are quantitatively interpreted by a simple 
model that yields reasonable values of the tip and metal 
work functions and of the tip/sample distance. The tun- 
neling dark current is described by a standard process 
between the semiconductor valence band and the metal- 
lic density of states. In contrast, the tunneling photocur- 
rent behavior is explained by a electric-field-dependent 
FN-like process between the semiconductor conduction 
band and the metallic density of states, including im- 
age charge effects. This mechanism accounts for all the 
results obtained before mechanical tip/sample contact : 
the exp(l/V) dependence of the tunnel photocurrent at 
large tip/sample distance, as well as the excess current at 
large tip bias observed just before contact. Once in con- 
tact the tunnel barrier for photoelectrons vanishes, and 
current flow is determined by the probablility of ballistic 
emission over the barrier. This results in a photocurrent 
dependent exponentially on V . 



For spin injection applications with GaAs tip: 



■5JL7 



it is 

possible that at low bias, injection will be best described 
by the standard-tunneling-based Nijmegen model. How- 
ever, the present work indicates that FN-like tunneling 
may also be observed at higher tip bias, even with GaAs. 
FN-like tunneling is likely to be spin dependent since the 
vacuum states into which photoelectrons tunnel should 
be hybridized with the spin polarized states of the mag- 
netic surface. In this case, FN-like tunneling would then 
be the relevant process (for example) for studies of the 
highly-polarized minority spin d-states of Fe and Co sur- 
faces which lie about 2 eV above the Fermi level 47 . 

In common with many works on field emission from 
silicon tips and on transport through Si/SiC>2 structures, 
current instabilities are observed. However, since the in- 
stabilities tend to occur over timescales longer than that 
taken for the spectral measurements of both the tunnel- 
ing dark current and the tunneling photocurrent, a mean- 
ingful analysis of the spectra is possible. With the aid 
of the model developed here, it is possible to show that 
the current instabilities are mostly related to changes in 
the thickness of the oxide layer covering the tip apex. 
Changes in this thickness may be related to electric field 
induced adsorption or desorption of foreign species such 
as water or oxygen. For this reason, an appropriate sur- 
face passivation of the tip is essential for future spin- 
polarized injection investigations using GaAs tips in liq- 
uid or neutral gas environments. Studies are underway 
to optimize such a tip treatement^. 
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APPENDIX A: NEGLIGIBLE FIELD 
ENHANCEMENT AT THE TIP APEX 

When discussing field emission phenomena, it is usual 
to describe the effect of tip geometry on the electric 
field by introducing the enhancement factor 7 defined by 
F ap ex = 7 -j j where F apex is the electric field at the tip 
apex. There is currently great interest in the use of sharp 



tips and high aspect ratio objects such as carbon nan- 
otubes, where field enhancement factors up to 7 = 10 4 
are reported^, since these make for very efficient field 
emitters. In the present case 7 is of order unity despite 
the fact that sharp tips are used. The main difference 
here is that the tip/sample distance d is of order 1 run 
(i.e. smaller than the tip radius, r sa 60 nm), whereas 
in field emission studies the cathode/anode distance can 
be as large as several /zm (i.e. much larger than the tip 
radius). At the scale of 1 nm therefore, the "sharp" tips 
used in this study are locally flat and 7 f=a 1. 

The electric field F(z) has been calculated along the 
axis between the tip apex and the surface. This calcula- 
tion was performed numerically in two dimensions, using 
a finite element resolution of Laplace's equation, with a 
fixed potential difference V between tip and sample. The 
tip is modelled as a conical section of half angle 12°, sim- 
ilar to the value for the tips used in the experiments, 
and the tip apex is terminated by a circle of radius r. 
As shown in Fig. [11] for a — 10 3 , which corresponds to 
a typical field emission experiment, the electric field is 
greatly enhanced near the tip surface where 7 sa 8. For 
a = 10 which corresponds to the present tunneling ex- 
periments, F(z) is constant and 7 = 1. The variation in 
7 with a is shown in Fig. Illf inset). It is noted that 7 
only deviates from unity for a larger than 1. 
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FIG. 11: Variation in normalized electric field (units of V/d) 
along the axis from the tip (z = 0) to the sample surface (z = 
1 reduced unit). A large increase in electric field is observed 
close to the tip for large a. The inset shows the variation 
in 7 with the aspect ratio, a = d/r . For small values of a, 
corresponding to the experiments reported here, 7 — » 1. 



